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Abstract 

In this article, we take the vector charmonium-like state y(4660) as a ip' fo{98Q) 
bound state (irrespective of the hadro-charmonium and the molecular state) tenta- 
tively, study its mass using the QCD sum rules, the numerical value My = 4.71 ± 
0.26 GeV is consistent with the experimental data. Considering the SU{3) symmetry 
of the light flavor quarks and the heavy quark symmetry, we also study the bound 
states V'o-(400 - 1200), T"7o(980) and T"'cr(400 - 1200) whh the QCD sum rules, 
and make reasonable predictions for their masses. 

PACS number: 12.39.Mk, 12.38.Lg 
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1 Introduction 

Two resonant structures are observed in the TT'^n^ip' invariant mass distribution in the 
cross section for the process e~^e~ — > Tr+Tr"'^' between threshold and ^/s = 5.5 GeV 
using 673 fb^^ of data on and off the T(45') (T'") resonance collected with the Belle 
detector at KEK-B, one at 4361 ± 9 ± 9 MeV with a width of 74 ± 15 ± 10 MeV, and 
another at 4664 ± 11 ± 5 MeV with a width of 48 ± 15 ± 3 MeV (they are denoted as 
y(4360) and y(4660) respectively), where the mass spectrum is parameterized by two 
Breit-Wigner functions [1]. The structure y(4660) is neither observed in the initial state 
radiation (ISR) process e+e" 7/s_R'?r"'"7r~ J/V' [2], nor in the exclusive cross processes 
e+e" ^ DD,DD*,D*D*,DDtt, J/i^D^*W^*^ [S H [SI El [7] . 

There have been several canonical charmonium interpretations for the y(4660), such 
as the S^^i state [SI, the 6^ Si state [9], the 5^ Si - A^Di mixing state [lO]. In Ref.[ll], 
Qiao suggests that the y(4660) is a baryonium state, the radial excited state of the 
^(|AcAc) + |ScS°)). In Ref.[I2j, Albuquerque et al take the y(4660) as a vector cscs 
tetraquark state, and study its mass with the QCD sum rules. A critical information 
for understanding the structure of those charmonium-like states is wether or not the vrvr 
comes from a resonance. There is some indication that only the y(4660) has a well 
defined intermediate state which is consistent with the scalar meson /o(980) in the mr 
invariant mass spectra [13j. In Ref.[14j. Guo et al take the y(4660) as a ^'/o(980) bound 
state (molecular state) considering the nominal threshold of the ^p' — /o(980) system is 
about 4666 ± 10 MeV [15], the y(4660) decays dominantly via the decay of the scalar 
meson /o(980), y(4660) V7o(980) iI^'ttti, ip'KK, the difficulties in the canonical 
charmonium interpretation can be overcome. In Refs.[T6l I17j . Voloshin et al argue that 
the charmonium-like states y(4660), Z(4430), y(4260), • • • may be hadro-charmonia. 
The relatively compact charmonium states {J/ip, ip' and Xcj) can be bound inside light 
hadronic matter, in particular inside higher resonances made from light quarks and (or) 
gluons. The charmonium state in such binding retains its properties essentially, the bound 
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system (hadro-charmonium, a special molecular state) decays into light mesons and the 
particular charmonium. 

In this article, we study the mass of the y(4660) as a ip' fo{980) bound state (irrespec- 
tive of the hadro-charmonium and the molecular state) using the QCD sum rules [18^ I19j. 
As a byproduct, we take into account the SU (3) symmetry of the light flavor quarks and 
the heavy quark symmetry, study the related hidden charm and hidden bottom states. In 
the QCD sum rules, the operator product expansion is used to expand the time-ordered 
currents into a series of quark and gluon condensates which parameterize the long distance 
properties of the QCD vacuum. Based on the quark-hadron duality, we can obtain copious 
information about the hadronic parameters at the phenomenological side [181 119j . 

The article is arranged as follows: we derive the QCD sum rules for the vector 
charmonium-like state y(4660) and the related bound states in section 2; in section 3, 
numerical results and discussions; section 4 is reserved for conclusion. 

2 QCD sum rules for the y(4660) and related bound states 

In the following, we write down the two-point correlation functions n^,y(p) in the QCD 
sum rules. 



where the Q denotes the heavy quarks c and b. We use the currents J^{x) and ti^{x) 
{Q = c) to interpolate the bound states V'/o(980) and ip'a{400 — 1200), respectively. 
The y(4660) can be tentatively identified as the '(/''/o(980) bound state, while there lack 
experimental candidates to identify the ip'a{400 — 1200) bound state. Considering the 
heavy quark symmetry, there maybe exist some hidden bottom bound states, for example, 
T/o(980), T7o(980), T'7o(980), T"7o(980), Tct(400- 1200), TV(400- 1200), T'V(400- 
1200), T'"(j(400 — 1200), we study those possibilities with the currents J^(x) and 
Vfiix) (Q = b), and make predictions for their masses which are fundamental parameters 
in describing a hadron. 

The hidden charm current c{x)j^c{x) can interpolate the charmonia J/tp, ip' , ^(3770), 
■0(4040), ■(/'(4160), 0(4415), while the hidden bottom current h{x)^^h{x) can inter- 
polate the bottomonia T, T', T", T'", T"", • • • |15j . We assume that the scalar mesons 
/o(980) and ct(400 — 1200) are the conventional qq states, to be more precise, they have 
large qq components. The currents J^{x) and r}^{x) {Q = c) have non- vanishing cou- 
plings with the bound states 7/0/0(980), 07o(98O), V;'7o(980), • • • and J/0(j(4OO- 1200), 
V'V(400 - 1200), 0"cj(4OO - 1200), • • • , respectively. The colored objects (diquarks) in a 
confining potential can result in a copious spectrum, there maybe exist a series of orbital 
angular momentum excitations; while the colorless objects (mesons) bound by a short 
range potential (through meson-exchange) should have a very limited spectrum, it is rel- 
atively easy to identify the molecule type bound states. We determine the masses of the 
ground states by imposing the two criteria of the QCD sum rules, then compare them 




(1) 



Q{x)-f^Q{x)s{x)s{x) , 
-^Q{x)^^Q{x) [u{x)u{x) +d{x)d{x 



)] 



(2) 



2 



with the nominal thresholds of the corresponding systems J/il^ — /o(980), ip' — /o(980), 
• • • . In Ref . [TB] , Voloshin et al argue that a formation of hadro-charmonium is favored for 
higher charmonium resonances ip' and XcJ as compared to the lowest states J and r/c. 

We can insert a complete set of intermediate hadronic states with the same quantum 
numbers as the current operators J^{x) and rj^{x) into the correlation functions Ii^y{p) to 
obtain the hadronic representation [18], |T9] . After isolating the ground state contributions 
from the pole terms of the Y and Z, we get the following result, 
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where the pole residues (or coupling) Ay and \z are defined by 

Xye^ = (0|J/7?^(0)|y(p)), 

AzPm = (0|J/r?^(0)|Z(p)), (4) 

and the is the polarization vector. In Eq.(3), we show the contribution from the scalar 
bound state Z explicitly, because the vector currents J^{x) and r/^(x) are by no means 
conserved. 

After performing the standard procedure of the QCD sum rules, we obtain two sum 
rules for the cess and hbss channels respectively (In the isospin limit, the interpolating 
currents result in two distinct expressions for the correlation functions n^jy(p), which are 
characterized by the number of the s quark they contain, thereafter will use the quark 
constituents to denote the corresponding quantities.): 

Aye m2" = / dsp{s)e , (5) 
J A 

r 1 GG 

pis) = Po(s) +P(ss>(s) + P{GG)i^) + P{GG)i^) i^:^) + P{ss)2is) . (6) 

The explicit expressions of the spectral densities Pois), P(ss)(s), pfGG)i^)' PfGG)i^) 
P{ss)^is) are presented in the appendix. The sq is the continuum threshold parameter 

„ l+Jl-4ml/s l-Jl-4m1/s --2 



and the M is the Borel parameter; a/ = — - — § ' ~ — — 2 ' 
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for the ccqq and bbqq channels with a simple replacement nig — > 0, (ss) — > (qq) and 
(sgsCrGs) {qgsoGq). 

We carry out the operator product expansion (OPE) to the vacuum condensates adding 
up to dimension-10. In calculation, we take assumption of vacuum saturation for high di- 
mension vacuum condensates, they are always factorized to lower condensates with vacuum 
saturation in the QCD sum rules, factorization works well in the large N^. limit. More- 
over, we neglect the terms proportional to the m„ and m^, their contributions are of minor 
importance. 

Differentiate the Eq.(5) with respect to then eliminate the pole residue Ay, we 
can obtain the sum rule for the mass of the bound state Y , 

Ja dspis)e ^ 
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3 Numerical results and discussions 



The input parameters are taken to be the standard values (qq) = —(0.24 it 0.01 GeV)^, 
{ss) = (0.8±0.2)(gg), {qgscrGq) = ml{qq), {sgsaGs) = ml{ss), = (0.8 ± 0.2) GeV^, 
(°^) = (0.33 GeV)^, rus = (0.14 ± 0.01) GeV, mc = (1.35 ± 0.10) GeV and = 
(4.8 ± 0.1) GeV at the energy scale /u = 1 GeV [H [M [20]. 

In the conventional QCD sum rules [18^ I19j . there are two criteria (pole dominance 
and convergence of the operator product expansion) for choosing the Borel parameter 
and threshold parameter sq. We impose the two criteria on the charmonium-like states 
Y to choose the Borel parameter and threshold parameter sq. The light tetraquark 
states cannot satisfy the two criteria, although it is not an indication non-existence of 
the light tetraquark states (For detailed discussions about this subject, one can consult 
Refs.[2ll|22]). 

We take the vector charmonium-like state y(4660) as the 'i/''/o(980) bound state ten- 
tatively, and take the threshold parameter as s^g = (4.66 -|- 0.5)^ GeV^ ?s 27GeV^ to take 
into account possible contribution from the ground state, where we choose the energy gap 
between the ground state and the first radial excited state to be 0.5 GeV. Taking into ac- 
count the SU{3) symmetry of the light flavor quarks, we expect the threshold parameter 
Sgg (for the bound state ^p'a{AOO — 1200)) is slightly smaller than the s%. Furthermore, 
we take into account the mass difference between the c and b quarks, the threshold pa- 
rameters in the hidden bottom channels are tentatively taken as s^^ = 144 GeV^ and 
s% = 145 GeV^. 

In this article, we take it for granted that the energy gap between the ground state and 
the first radial excited state is about 0.5 GeV, and use this value as a guide to determine 
the threshold parameter sq with the QCD sum rules. 

The contributions from the high dimension vacuum condensates in the operator prod- 
uct expansion are shown in Figs. 1-2, where (and thereafter) we use the {qq) to denote the 
quark condensates {qq), {ss) and the {qgsCrGq) to denote the mixed condensates {qgsCrGq), 
{sgsCrGs). From the figures, we can see that the contributions from the high dimension con- 
densates change quickly with variation of the Borel parameter at the values < 2.8 GeV^ 
and < 7.5 GeV^ for the hidden charm and hidden bottom channels respectively, such 
an unstable behavior cannot lead to stable sum rules, our numerical results confirm this 
conjecture, see Fig. 4. 

At the values > 2.8 GeV^ and sq > 26GeV^, the contributions from the {qq)'^ + 
{qq){qgsCrGq) term are less than (or equal) 18.5% for the cess channel, the corresponding 
contributions are less than (or equal) 36.5% for the ccqq channel; the contributions from the 
vacuum condensate of the highest dimension {qgscGq)'^ are less than 5% for all the hidden 
charm channels, we expect the operator product expansion is convergent in the hidden 
charm channels. At the values > 7.6 GeV^ (In Figs. 2-4, the vertical line corresponds 
to the value = 7.6 GeV^ in the hidden bottom channels.) and sq > 148 GeV^, the 
contributions from the {qq)'^ + {Q'l){Q9s(^Gq) term are less than 7% for the bbss channel, 
the corresponding contributions are less than (or equal) 18% for the bbqq channel; the 
contributions from the vacuum condensate of the highest dimension {qggCrGq)'^ are less 
than (or equal) 7% for all the hidden bottom channels, we expect the operator product 
expansion is convergent in the hidden bottom channels. 

The contributions from the gluon condensate ( °''^'^ ) are rather large, while the contri- 
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butions from the high dimension condensates (^^^^) [{qq) + {qgscGq) + (qq)'^] are small 
enough, the total contributions involving the gluon condensate are less than (or equal) 
30% (22%) for the cess (ccqq) channel at the values > 2.8 GeV^ and sq > 26GeV2; 
while the contributions are less than 21% (17%) for the bbss (bbqq) channel at the values 
M2 > 7.6 GeV^ and sq > 148 GeV^. In the QCD sum rules for the tetraquark states 
(irrespective of the molecule type and the diquark-antidiquark type), the contributions 
from the gluon condensate are suppressed by large denominators and would not play any 
significant roles for the light tetraquark states |231 124j . the heavy tetraquark state [21] 
and the heavy molecular state ^25j; the present sum rules seem rather exotic. If we take 
a simple replacement s{x)s{x) — > {ss) and [u{x)u{x) + d{x)d{x)~\ — > 2{qq) in the interpo- 
lating currents J^{x) and ^/^(x), we can obtain the standard vector heavy quark current 
Q{x)'yfj,Q{x), where the gluon condensate plays an important rule in the QCD sum rules 

m- 

In calculation, we observe that the dominant contributions come from the perturbative 
term and the (qq) + {qggCrGq) term at the values > 2.8 GeV^ and sq > 26 GeV^ for 
the hidden charm channels and at the values > 7.6 GeV^ and sq > 148 GeV^ for the 
hidden bottom channels, the operator product expansion is convergent. 

In this article, we take the uniform Borel parameter M^^„, i.e. M^-^ > 2.8 GeV^ and 
^min — '^■^ GeV^ for the hidden charm and hidden bottom channels, respectively. 

In Fig. 3, we show the contributions from the pole terms with variation of the Borel pa- 
rameters and the threshold parameters sq- If the pole dominance criterion is satisfied, 
the threshold parameter sq increases with the Borel parameter monotonously. From 
Fig.3-A, we can see that the pole dominance criterion cannot be satisfied at the values 
So < 25 GeV^ and > 2.8 GeV^ in the ccss channel, the threshold parameter sq has to 
be pushed to larger value. 

The pole contributions are larger than 45% at the values < 3.2 GeV^ and sq > 
25GeV2, 26GeV2 for the ccqq, ccss channels respectively; and larger than 50% at the 
values < 8.2GeV^ sq > 146 GeV^ 148 GeV^ for the bbqq and bbss channels respec- 
tively. Again we take the uniform Borel parameter M^^^, i.e. M^^^ < 3.2 GeV^ and 
^max — GeV^ for the hidden charm and hidden bottom channels, respectively. 

If we take uniform pole contributions, the interpolating current with more s quarks 
requires larger threshold parameter due to the SU{3) breaking effects, see Fig. 3. The 
threshold parameters in the ccqq and bbqq channels are slightly smaller than the ones in the 
ccss and bbss channels respectively. In this article, the threshold parameters are taken as 
So = (_26±1) GeV^ (27±1) GeV^, (148±2) GeV^ and (150±2) GeV^ for the ccqq, ccss, bbqq 
and bbss channels, respectively; the Borel parameters are taken as = (2.8 — 3.2) GeV^ 
and (7.6 — 8.2) GeV^ for the hidden charm and hidden bottom channels, respectively. 
In those regions, the pole contributions are about (45 — 69)%, (46 — 69)%, (50 — 66)% 
and (51 — 67)% for the ccss, ccqq, bbss and bbqq channels, respectively; the two criteria 
of the QCD sum rules are fully satisfied [18l [19]. Naively, we expect the bound state 
with the scalar meson /o(980) will have larger mass than the corresponding one with the 
scalar meson (j(400 — 1200), our numerical calculations confirm this conjecture, see Fig. 4. 
Although smaller threshold parameters lead to slower convergent behavior in the operator 
product expansion, the two criteria of the QCD sum rules are still satisfied, one can consult 
Figs. 1-2. 

The Borel windows Af^^^, — M^^^ change with variations of the threshold parameters 
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Figure 1: The contributions from different terms with variation of the Borel parameter 
in the operator product expansion. The A, B, C and D correspond to the contribu- 
tions from the (qgsCrGq)'^ term, the (qq)'^ +{QQ){Q9sO'Gq) term, the ( ) term and the 
(^^) + (^^) [{qq) + (qgsCrGq) + {qqf] term, respectively. The (I) and (II) denote the 
cess and ccqq channels, respectively. The notations a, (3, 7, A, p and r correspond to the 
threshold parameters sq = 23GeV2, 24GeV2, 25GeV2, 26GeV2, 27GeV2 and 28GeV2, 
respectively. 
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Figure 2: The contributions from different terms with variation of the Borel parameter 
in the operator product expansion. The A, B, C and D correspond to the contribu- 
tions from the (qgsCrGq)'^ term, the (qq)'^ +{QQ){Q9sO'Gq) term, the ( ) term and the 
(^^) + (^^) [{qq) + (qgsCrGq) + {qqf] term, respectively. The (I) and (II) denote the 
bbss and bbqq channels, respectively. The notations a, /3, 7, A, p and r correspond to 
the threshold parameters sq = 142GeV^ 144 GeV^ 146GeV^ 148 GeV^ 150 GeV^ and 
152 GeV^, respectively. 



7 





Figure 3: The contributions from the pole terms with variation of the Borel parameter 
M2. The A, B, C, and D denote the cess, ccqq, bbss and bbqq channels, respectively. In 
the hidden charm channels, the notations a, /3, 7, A, p and r correspond to the threshold 
parameters sq = 23GeV2, 24GeV^ 25 GeV^, 26Ge\^, 27Ge\'^ and 28GeV^ respectively 
; while in the hidden bottom channels they correspond to the threshold parameters sq = 
142 GeV^, 144 GeV^, 146 GeV^, 148 GeV^, 150 GeV^ and 152 GeV^, respectively. 
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Figure 4: The masses of the vector bound states with variation of the Borel parameter 
M2. The A, B, C, and D denote the cess, ccqq, bbss, and bbqq channels, respectively. In 
the hidden charm channels, the notations a, /3, 7, A, p and r correspond to the threshold 
parameters sq = 23GeY'^, 24GeV2, 25GeY^, 26GeV^ 27GeV2 and 28GeV^ respectively 
; while in the hidden bottom channels they correspond to the threshold parameters sq = 
142 GeV^, 144 GeV^, 146 GeV^, 148 GeV^, 150 GeV^ and 152 GeV^ respectively. The ^ 
and fi denote the J/'0— /o(980) and ip' — fo{980) thresholds respectively in the cess channel, 
while in the bbss channel they correspond to T" — /o(980) and T'" — /o(980) thresholds 
respectively. 
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Figure 5: The masses of the vector bound states with variation of the Borel parameter 
M^. The A, B, C, and D denote the cess, ccqq, bbss, and bbqq channels, respectively. 
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Figure 6: The pole residues of the vector bound states with variation of the Borel param- 
eter M^. The A, B, C, and D denote the cess, ccqq, bbss, and bbqq channels, respectively. 
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So, see Fig. 3. In this article, the Borel windows are taken as 0.4 GeV^ and 0.6 GeV^ 
for the hidden charm and hidden bottom channels respectively, they are small enough. 
Furthermore, we take uniform Borel windows and smear the dependence on the threshold 
parameters so- If we take larger threshold parameters, the Borel windows are larger and 
the resulting masses are larger, see Fig. 4. In this article, we intend calculate the possibly 
lowest masses which are supposed to be the ground state masses by imposing the two 
criteria of the QCD sum rules. 

In Fig. 4, we plot the bound state masses My with variation of the Borel parame- 
ters and the threshold parameters. The hidden charm current c(x)7^c(x) can interpo- 
late the charmonia J/ip, ip' , V'lSTTO), 7/;(4040), ^^(4160), ^^(4415), while the hid- 
den bottom current b{x)j^b{x) can interpolate the bottomonia T, T', T", T'", T"", • • • 
|15j . The currents have non- vanishing couplings with the bound states J/'i/'/o(980), 

V'7o(980), V'7o(980), ••• and T/o(980), T7o(980), T'7o(980), T"7o(980), respec- 
tively. From Figs.3-A,3-C,4-A,4-C, we can see that the QCD sum rules support existence 
of the V7o(980) and T"7o(980) bound states, the nominal thresholds of the J/ip- fo{980) 
and T" — /o(980) systems are too low, and we cannot reproduce the J /^pf o{980) and 
'^'7o(980) bound states. Our numerical results support the conjecture of Voloshin et al, 
i.e. a formation of hadro-charmonium is favored for higher charmonium resonances ip' and 
Xcj as compared to the lowest states J/tp and rjc [16] . 

In this article, we intend prove that the -070(980) and T"7o(980) bound states can be 
reproduced by the QCD sum rules, the charmonium-like state y(4660) has the possibility 
to be a ip' fo{980) bound state. 

Taking into account all uncertainties of the input parameters, finally we obtain the 
values of the masses and pole resides of the vector bound states Y, which are shown in 
Figs. 5-6 and Tables 1-2. In this article, we calculate the uncertainties 6 with the formula 



i ^ * ^ 

where the / denote the hadron mass My and the pole residue Ay, the Xi denote the input 
QCD parameters ttt-c, nib, {qq), (ss), • • • . As the partial derivatives ^ are difficult to carry 

out analytically, we take the approximation (^§i'^ (xj — Xj)^ [fixi ± Axj) — /(x,)]^ in 
the numerical calculations. 

From Tables 1-2, we can see that the uncertainties of the masses My are rather small 
(about 5% in the hidden charm channels and 2% in the hidden bottom channels), while the 
uncertainties of the pole residues Ay are rather large (about (30 — 50)%). The uncertainties 
of the input parameters {{qq), (ss), {sgsCrGs), {qgsCrGq), nis, rric and m^) vary in the 
range (2 — 25)%, the uncertainties of the pole residues Ay are reasonable. We obtain the 
squared masses My through a fraction, see Eq.(7), the uncertainties in the numerator 
and denominator which origin from a given input parameter (for example, (ss), {sgsCrGs)) 
cancel out with each other, and result in small net uncertainty. 

In table 1, we also present the nominal thresholds of the V'' — /o(980), tp' — a{400— 1200), 
T'" - /o(980) and V" - cj(400 - 1200) systems. From the table, we can see that the 
y(4660) can be tentatively identified as the "070(980) bound state. The predicted mass of 
the •0'o-(4OO - 1200) bound state is about (4.59 ± 0.19)GeV, while the nominal threshold 
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bound states 


My (GeV) 


M^,/T"' + ^f„/<x (GeV) 


My (GeV)* 


cess 


4.71 ±0.26 


4.666 


4.63 


ccqq 


4.59 ±0.19 


4.086 - 4.886 


4.56 


bbss 


11.57 ±0.20 


11.559 


11.56 


bbqq 


11.52 ±0.18 


10.979 - 11.779 


11.51 



Table 1: The masses of the bound states, we use the star * to denote the central values 
from the sum rules where the perturbative contributions are multiplied by a factor 2. 



bound states 


Ay (10-^GeV^) 


Ay (10-^GeV^)* 


cess 


3.70+}-^« 


5.23 


ccqq 


3.49li:- 


4.84 


bbss 


19.2 ±8.2 


28.6 


bbqq 


19.2 ±6.7 


27.3 



Table 2: The pole residues of the bound states, we use the star * to denote the central 
values from the sum rules where the perturbative contributions are multiplied by a factor 
2. 

of the ip' - cr(400 - 1200) system is about (4.086 - 4.886) GeV. There maybe exist such a 
bound state. The V'''?" (400— 1200) bound state can be produced in the initial state radiation 
process e^e~ 'JiSR'^^'^~'4'' or in the exclusive decays of the B meson through b ccq 
at the quark level. There still lack experimental candidates to identify the ■0'o"(4OO — 1200) 
bound state, such a bound state is difficult to observe due to the broad width of the scalar 
meson cj(400 - 1200). 

In the bbss channel, the numerical result My = 11.57 ± 0.20 GeV indicates that there 
maybe exist a T"'/o(980) bound state, which is consistent with the nominal threshold 
Mr'" + Mfg = 11.559 GeV, while the nominal thresholds Mr + Mj^ = 10.44 GeV, Mx' ± 
Mf^ = 11.00 GeV, Mrii+Mf^ = 11.335 GeV are too low. The scalar meson cr(400- 1200) 
is rather broad with the Breit-Wigner mass formula (400 — 1200) — i(250 — 500) |15j . 
Considering the SU{3) symmetry of the light flavor quarks, we can obtain the conclusion 
tentatively that there maybe exist the V''o"(400 - 1200) and T"'cr(400 - 1200) bound states 
which lie in the regions (4.086 - 4.886) GeV and (10.979 - 11.779) GeV, respectively. As 
the energy gaps between the T's are rather small and the scalar meson (t(400 — 1200) is 
broad enough, there maybe exist the Tcr (400- 1200), T'c7(400-1200) and T'V (400- 1200) 
bound states. We cannot draw decisive conclusion with the QCD sum rules alone. 

At the energy scale /i = IGeV, ^ 0.19 [26j, if the perturbative 0{as) corrections 
to the perturbative term are companied with large numerical factors, 1 ± ^(s,rn-Q)^, for 
example, ^(s, mg) > ^ ~ 5, the contributions may be large. We can make a crude estima- 
tion by multiplying the perturbative term with a numerical factor, say 1 ±^(s, mg)^ = 2, 
the masses My decrease slightly while the pole residues Ay increase remarkably, see Ta- 
bles 1-2. The main contribution comes from the perturbative term, the large corrections 
in the numerator and denominator cancel out with each other (see Eq.(7)). In fact, the 
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^(s,mQ) are complicated functions of the energy s and the mass mg, such a crude esti- 
mation maybe underestimate the 0{as) corrections, the uncertainties originate from the 
0{as) corrections maybe larger. 

The charmonia J/ip, ip', -0(3770), V(4040), -0(4160), -0(4415), • • • and the bottomonia 
T, T', T", T'", T"", • • • also have Fock states with additional qq components beside the 
QQ components. The currents Jfi{x) and 'r]^(x) may have non- vanishing couplings with 
the charmonia and bottomonia, those couplings are supposed to be small, as the main 
Fock states of the charmonia and bottomonia are the QQ components, and the charmonia 
and bottomonia have much smaller masses than the corresponding molecular states Y. 

In this article, we take the assumption that the scalar mesons /o(980) and fT(400— 1200) 
are the conventional qq mesons, or more precise, they have large qq components. There 
are hot controversies about their nature, for example, the conventional qq states (strongly 
affected by the nearby thresholds), the tetraquark states, the molecular states \27\ I28j. In 
Ref . [29] ■ we take the scalar mesons ao(980) and /o(980) as the conventional qq mesons, 
study the strong couplings to the nearby thresholds, and observe that the strong couplings 
are rather large. Then we draw the conclusion that the ao(980) and /o(980) may have a 
small qq kernel of the typical qq meson size, strong coupling to the nearby KK threshold 
may result in some tetraquark components (irrespective of a nucleon-like bound state and 
a deuteron-like bound state) [29]. The decay /o(980)/cr(400 — 1200) — > vt-tt, KK can occur 
through the tetraquark quark components naturally. 

The LHCb is a dedicated b and c-physics precision experiment at the LHC (large 
hadron collider). The LHC will be the world's most copious source of the b hadrons, 
and a complete spectrum of the b hadrons will be available through gluon fusion. In 
proton-proton collisions at ^/s = 14TeV, the bb cross section is expected to be ~ bOOfib 
producing 10^^ bb pairs in a standard year of running at the LHCb operational luminosity 
of 2 X 1032cjj^-2ggp-i [3Qj_ rp^g bound states T"7o(980) and T"V(400 - 1200) predicted 
in the present work may be observed at the LHCb, if they exist indeed. We can search for 
those bound states in the Ttttt, T'tttt, V'tttt, V'tttt, TKK, VKK, VKK, V'KK, ■ ■ ■ 
invariant mass distributions. 

4 Conclusion 

In this article, we take the the vector charmonium-like state y(4660) as the -0'/o(98O) 
bound state (irrespective of the hadro-charmonium and the molecular state) tentatively, 
study its mass using the QCD sum rules, the numerical result My = 4.71 it 0.26 GeV 
is consistent with the experimental data 4664 it 11 it 5 MeV. Considering the SU{3) 
symmetry of the light flavor quarks and the heavy quark symmetry, we also study the 
bound states -0V(4OO - 1200), T"7o(980) and T"V(400 - 1200) with the QCD sum rules, 
and make reasonable predictions for their masses. Our predictions depend heavily on the 
two criteria (pole dominance and convergence of the operator product expansion) of the 
QCD sum rules. We can search for those bound states at the LHCb, the KEK-B or the 
Fermi-lab Tevatron. 
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Appendix 

The spectral densities at the level of the quark-gluon degrees of freedom: 



3 Z*^/ f'l—a 

Pois) = 4Q9g^6 J y dl3aP{l- a- pf{s-m'^Qf{5s-fn^) 



10247r6 



da 



l-a 



df3{l-a- f3f{s-fn, 



2 \3 



P{ss){s) = 



1287r4 



da I dPaP{s — mQ){3s — mq) 
hi 



647r4 



uisisgsCrGs) f"f 



da I df3{s — nin) 
Pi 



327r4 



~2 

/ daa{l — a){2s — mq) 

Joii 



2 'msmQ{sgs(TGs) r^f 



327r4 
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J Oii 



da , 



nf ~2 rn^iss)'^ f^f 
P(ss)^is) = / daa{l - a){2s - niQ) -„ n / da 



327r2 
32^^2 
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3 + ( 3s + J S{s - ml) 



da 



1 + 



M2 



1287r2M2 
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1 + — + 



M2 2M4 



(5(s - mg) 



3m.Q{sgsaGs)^ f^^l ~2 
H — ? — / das 5(s — mn) , 
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3 f^f pi—ot 
P(GG>(«) = 1024^ / / dPaP{s - ml){3s - -ml) 

J Oti J 
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Jai J Pi 
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1 1 



m 



Q 



10247r4 

ms{ss) C'^f 
' 1287r2 

msmQ^ss) f^f 
' 1287r2M2 



1 1 



aj- fl~a 

da I dp 

Jpi 

t-l—a 

da I dp[Z + sS{s - rhn)] 

J(3r 

paf rl — a 

da dp 
Jai Jg 



{l-a-pf{2s-m^Q) 
{l-a-P)^{s-ml) 
(l-a-pf 



I Hi 



a P 
P^'^7^ 



msmQ{ss) /■"/ 



1287r2M2 



I da I dp 

Jai J Pi 



1287r2 



fl—a 

da I dp 

ai J Pi 



1 1 

1 1 ■ 



- friQ) 



(12) 



5767r2 



af 



da 



2 + 



M2 



288 



•Jai 



da 



2 + 



M2 



msW?Q{sgsaGs) r'^f 



5767r2M2 



/ 



da 



1 — a 



+ 



a 



a^ 



msUiQisgsCrGs) /""/ 
H / aa 



5767r2M4 
" 288M2 / 



1 



+ 



(1 - a)2 
1 



M2 



da 



+ ^ ' / da 
288M4 



1 — a 
1 



a^ (1 — a)3 

a 



+ 



(l-a)2 
1 



1 - 



M2 



a^ (1 — a)3 



rUsm'^isgsCrGs) f^f 
1927r2Af2 

96M2 



da 



da 



1 



1 



1 



a2 (1 — a)2 



a2 (l-a)2_ 



5{s - mq) . 



(13) 



16 



Acknowledgements 



This work is supported by National Natural Science Foundation of China, Grant Number 
10775051, and Program for New Century Excellent Talents in University, Grant Number 
NCET-07-0282. 

References 

[1] X. L. Wang et al, Phys. Rev. Lett. 99, 142002 (2007). 

[2] C. Z. Yuan et al, Phys. Rev. Lett. 99, 182004 (2007). 

[3] G. Pakhlova et al, Phys. Rev. Lett. 98, 092001 (2007). 

[4] G. Pakhlova et al, Phys. Rev. Lett. 100, 062001 (2008). 

[5] G. Pakhlova et al, Phys. Rev. D 77, 011103 (2008). 

[6] B. Aubert et al, larXiv:0710. 13711 

[7] P. Pakhlov et al, Phys. Rev. Lett. 100, 202001 (2008). 

[8] G. J. Ding, J. J. Zhu and M. L. Yan, Phys. Rev. D 77, 014033 (2008). 

[9] B. Q. Li and K. T. Chao. [a?Xw:0903.550 6' 

[10] A. M. Badahan, B. L. G. Bakker and I. V. Danilkin, larXiv:0805.2291l 

[11] C. F. Qiao, J. Phys. G35, 075008 (2008). 

[12] R. M. Albuquerque and M. Nielsen, Nucl. Phys. A 815, 53 (2009). 

[13] R. Faccini, larXiv:0801.2679l 

[14] F. K. Guo, C. Hanhart and Ulf-G. Meissner, Phys. Lett. B 665, 26 (2008). 

[15] C. Amsler et al, Phys. Lett. B667, 1 (2008). 

[16] S. Dubynskiy and M. B. Voloshin, Phys. Lett. B 666, 344 (2008). 

[17] M. B. Voloshin, Prog. Part. Nucl. Phys. 61 (2008) 455. 

[18] M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B147 (1979) 385, 
448. 

[19] L. J. Reinders, H. Rubinstein and S. Yazaki, Phys. Rept. 127 (1985) 1. 

[20] B. L. loffe. Prog. Part. Nucl. Phys. 56 (2006) 232. 

[21] Z. G. Wang, Eur. Phys. J. C62 (2009) 375. 

[22] Z. G. Wang, Chin. Phys. C32 (2008) 797. 

[23] Z. G. Wang, Nucl. Phys. A791 (2007) 106. 

17 



[24] Z. G. Wang, W. M. Yang and S. L. Wan, J. Phys. G31 (2005) 971. 

[25] Z. G. Wang, Eur. Phys. J. C63 (2009) 115. 

[26] M. Davier, A. Hocker, Z. Zhang, Rev. Mod. Phys. 78 (2006) 1043. 

[27] F. E. Close and N. A. Tornqvist, J. Phys. G28 (2002) R249. 

[28] C. Amslcr and N. A. Tornqvist, Phys. Kept. 389 (2004) 61. 

[29] Z. G. Wang, W. M. Yang and S. L. Wan, Eur. Phys. J. C 37, 223 (2004). 

[30] G. Kane and A. Pierce, "Perspectives On LHC Physics", World Scientific Publishing 
Company, 2008. 



18 



